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In part I of this paper a modified BCS mechanism of Cooper pair formation of electrons was 
proposed. This mechanism is connected with the existence of a narrow, roughly half-filled "super- 
conducting energy band" of given symmetry. The special symmetry of the superconducting band 
was interpreted within a nonadiabatic extension of the Heisenberg model of magnetism. Within 
this nonadiabatic Heisenberg model, the electrons of the superconducting band are constrained to 
form Cooper pairs at zero temperature because in any normal conducting state the conservation 
of crystal-spin angular momentum would be violated. Except for this participation of the angu- 
lar momentum, the pair formation is mediated in the familiar way by phonons. Superconducting 
bands can be identified even within a free-electron band structure. Therefore, in this paper the band 
structures of the bcc and hep solid solution alloys composed of transition elements are approximated 
by appropriate free-electron band structures with s-d symmetry. From the free-electron bands, the 
number n of valence electrons per atom related to the maxima of the superconducting transition 
temperature T c in these solid solutions is calculated within the nonadiabatic Heisenberg model. The 
observed two maxima of T c are reproduced without any adjustable parameter at valence numbers 
n approximately equal to 4.9 and 6.5 in bcc and 4.7 and 6.7 in hep solid solutions. This result is in 
good aggreement with the measured values of 4.7 and 6.4 of Hulm and Blaugher. The upper maxi- 
mum is predicted not to exist in bcc transition elements but to occur in several ordered structures 
of bcc solid solution alloys. 

PACS numbers: 74.20.-z, 74.25.Jb, 74.62.Bf, 74.70.Ad 



I. INTRODUCTION 



In an early paper, MatthiasS pointed out that the su- 
perconducting transition temperature T c depends cru- 
cially on the number n of valence electrons per atom. 
Matthias proposed that n — 5 and n — 7 valence elec- 
trons per atom are particularly favorable for supercon- 
ductivity. By means of this empirical rule, which is often 
referred to as Matthias rule, a great number of supercon- 
ducting alloys composed of transition as well as nontran- 
sition elements have been discovered. 

Later, Hulm and BlaugheiH presented experimental 
data on the variation of T c with composition over the en- 
tire range of binary body-centered cubic solid solutions 
in the transition metal alloy system. They found two 
distinct maxima of T c , one located between group 4 and 
group 5 close to the valence number n ~ 4.7, and the 
second located between group 6 and group 7 close to the 
valence number n s» 6.4; see Fig. |. Though the posi- 
tions of the observed maxima differ from the values of 
n = 5 and n — 7 originally proposed by Matthias, the re- 
sults of Hulm and Blaugher confirm the essential physical 
substance of the Matthias rule, namely that T c strongly 
depends on the number n of valence electrons per atom 
and that this n dependence is similar in every alloy sys- 
tem. 



12 ■ 



c 



6 ■ 



1 1 1 1 


— i — i — i — i — 




- Zr-Nb 




Ti V Cr Mn 










Zr-Nb-Mo^Tc 






Hexagonal J' 




Hf To W Re 




9 ' Sigma " 
/ 1 Region 








■Region / 




7 ! *' 


if Ti-Nb 




• / Ko«. 






V 


A Cubic 




-f Limit 




A Region 




■ / Limit 


A of 






1 of 


J Cubic 




/ Cubic 


* Phase 

J . . . 


1— — 1 1 1 


J Phase R« & 
i i I i 1 



u 

Ti 
Zr 



5 

Nb 



6 

Mo 



n 



7 

Re 

(Tc) 



FIG. 1. Transition temperature T c versus number n of va- 
lence electrons per atom in solid solution alloys of the transi- 
tion elements, as measured by Hulm and Blaugher (Ref. |^). 
The positions of the maxima are n « 4.7 and n « 6.4. 

As is well-known, there is overwhelming theoreti- 
cal and experimental evidence that an electron-phonon 
mechanism is responsible for superconductivity in the 
transition elements^ In particular, the variation of T c 
with the number n of valence electrons per atom is under- 
stood within the McMillan theoryo of.-superconductors 
with strong electron-phonon coupling lira Nevertheless, it 
cannot be excluded that this variation of T c can also be 
understood in a new way provided that this new way does 
not contradict the standard theory of superconductivity. 

In this paper, the Matthias rule is-jpiconsidered within 
the nonadiabatic Heisenberg modelQia which is defined 
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by three new postulates given, e.g., in Ref. | 

Niobium possesses a narrow, roughly half-filled energy, 
band of special symmetry which in part I of this papert2l 
was called "superconducting band" (cr band). It was 
shown that if the postulates of the nonadiabatic Heisen- 
berg model are satisfied within this a band, then the 
electrons of this band form Cooper pairs at zero temper- 
ature. pJust as in the Bardeen-Cooper-Schrieffer (BCS) 
theorylill, the coupling of the electrons is mediated by 
phonons. However, the formation of Cooper pairs is con- 
strained in an absolutely new way by the conservation 
law of spin angular momentum: within a narrow, partly 
filled a band the crystal-spin angular momentum is con- 
served at zero temperature if and only if the electrons 
form Cooper pairs. 

Within the nonadiabatic Heisenberg model, the super- 
conducting transition temperature T-! of the jth a band 
of a metal may be calculated in the BCS limit from the 
modified BCS equation 



(1.1) 



see Sec. IIIC of part I. V and 8 denote the effective 
electron-phonon interaction and the Debye temperature, 
respectively. The nonadiabatic Heisenberg model mod- 
ifies the interpretation of the BCS formula in that the 
parameter Nj(Ep) now is only a partial density of states 
of the electrons at the Fermi level, namely the density of 
states of the electrons of the jth a band. 

In part I of this paper, only bec metals with one atom 
per unit cell have been considered. The present paper, 
however, deals with materials with more than one atom 
in the unit cell. Therefore, single a bands must be re- 
placed by "cr-band complexes" containing just as many 
bands as the re ar e atoms in the unit cell; see appendix [B|. 
In equation (1.1) the index j now labels the jth cr-band 
complex. In this paper, j — L and j = U denote the 
cr-band complex belonging to the lower and upper maxi- 
mum of T c , respectively. 

In Sec. H free-electron energy bands will be con- 
structed which in this paper shall approximate the s-d 
energy bands of the binary bec and hep solid solutions 
in the transition metal alloy system. That means, a free- 
electron band structure is constructed consisting only of 
linear combinations of plain waves with the symmetry of 
a s-d band group. Because cr-band complexes can already 
be identified within these free-electron bands, in Sec. Ill 
the parameters N1(Ep) and N^(Ep) will be calculated 
within this free-electron model. With these iV£ u (Ep), 
the positions of the maxima of T C L and will be cal- 
culated by means of Eq. (1.1) without any fitting pa- 
rameter, whereas the constant V will be adjusted such 
that the maximum T c amounts to about 10 K . This ad- 
justment does not influence the positions of the maxima 
noticeably. 

In appendix |^ the definition of a cr-band complex will 
be given in a metal with any given space group and any 
given number of atoms per unit cell, and in appendix |y 



the related spin dependent Wannier functions (spin de- 
pendent Wfs) will be defined. 



II. FREE-ELECTRON ENERGY BANDS WITH 
s-d SYMMETRY 

The energy bands of the valence electrons of the tran- 
sition elements are often referred to as s-d bands since 
the valence electrons in isolated atoms of these metals 
have s or d symmetry. Consider a metal with the space 
group G and [l atoms in the unit cell. In this paper, 6/i 
energy bands of this metal are called "s-c? bands" if, first, 
the Bloch functions of these bands can be unitarily trans- 
formed into Wfs centered at the atoms and if, secondly, 
at each atom there are six Wfs with one of them having 
s-like and the other five having rf-like symmetry, s- and 
d-like means that in Eq. (1.8) of Ref. O the Wfs belong 
to a representation of G which is compatible with the 
representation Dq + Di of the three-dimensional rotation 
group 0(3). 

The band structure of the transition metals shall be 
approximated by a set of free-electron bands which sat- 
isfy the following three conditions. 

(1) The set shall consist of 5/Lt bands with the remaining 
[i bands connecting it with higher lying, bands. 

(2) Within the representation domainEJ of the Brillouin 
zone the energy Ei (k) of the ith band shall be calculated 
by 



Ei(k) = — — (k — Km) 2 , 
Am 



(2.1) 



where Km denotes the wave vector of the reciprocal lat- 
tice point M related according to Table | to the ith band. 
(3) The bands shall be labeled by the representations 
given in Table |n] and III, respectively. Since, however, 
the set consists of 5fi bands only, some representations 
may be absent. 



A. Body-centered cubic structure 

In the bec structure we have [i = l. Table gives all 
the representations of a set of six s-d bands in the Bril- 
louin zone of the bec lattice as derived from Table 2.7 of 
Ref. || 

In the band structure of the typical transition element 
Nb [the calculation of MattheialJ is depicted, e.g., in 
Fig. 1 of part I of this paper], the majority of the Bloch 
functions has the symmetry given in Table P However, 
Nb and the other transition elements do not possess an 
ideal set of s-d bands because one level near the Fermi 
energy at point N is labeled by the p-like representation 
N{. ' 

The symmetry of the free-electron wave functions 
(which are symmetrized plane waves) ca«ibe determined 
by standard group theoretical methods. 
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It turns out that no set of free-electron bands exists 
which complies with all the three conditions and which 
satisfies the compatibility relations in the Brillouin zone. 
Therefore, as suggested by the band structures of the 
bcc transition elements, I replace one of the two N± lev- 
els by the p-like N[ level. With this modification of the 
third condition we get the one free-electron band struc- 
ture given in Fig. |^ which is used in this paper as approx- 
imation for the energy bands of the bcc solid solution 
alloys of the transition elements. 
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FIG. 2. Free-electron s-d energy bands which serve as ap- 
proximation of the energy bands of the bcc solid solution 
alloys of the transition elements. The bands are numbered 
arbitrarily from 1 to 6 with the related points M of the re- 
ciprocal lattice being given in Table |a. Band 3 (dotted line) 
is a superconducting band. Band 4 (bold line) does not form 
a superconducting band, but it becomes a superconducting 
band complex when it is folded into the Brillouin zone of 
the MoSi2 structure (see Fig or into the Brillouin zone of 
several atomic arrangements with the hexagonal space group 
D 3 

The symmetry structure of the free-electron bands has 
a close resemblance to what is found in the band struc- 
tures of the bcc transition elements, whereas, of course, 
the free-electron energy Ei(k) deviates considerably from 
the actual energy. By "symmetry structure" I mean the 
occurence and the connections of the representations be- 
tween each other. For instance, in the free-electron case 
as well as in the actual bands, the H\2 state occures and 
is connected with the P4 and the iVi and N[ states. 



B. Hexagonal close-packed structure 



rotation group 0(3).] 

The free-electron band structure which can be con- 
structed by starting from the three conditions at the top 
of this section, is depicted in Fig. |[ Other sets of free- 
electron bands which also comply with these conditions, 
only differ unessentially from the band structure depicted 
in Fig. [3| For instance, the upper Tf and states may 
be interchanged. Hence, I use the band structure in Fig. y 
as approximation of the energy bands of the hep solid so- 
lution alloys of the transition elements. 

In the case of the hep structure, it is not necessary 
to use p-like states in the free-electron s-d band struc- 
ture in order for the compatibility relations to be ful- 
filed. In the real band structures of the hep transition 
elements, however, p-like states are present and influence 
the symmetry-structure of the bands. In rhenium, e.g., a 
p-like state, namely , is situated near the Fermi level; 
see Fig. 2 of Ref. [It]. Therefore, an additional band exists 
in the band structure which causes a change of the sym- 
metry structure. In this paper, this pertubation by the 
p-electrons will be ignored since it is difficult to include 
p-like states in an unequivocal way. 



III. CALCULATION OF THE POSITIONS OF 
THE MAXIMA OF T c . 

In this section I shall calculate the positions of the 
maxima of the su perc onducting transition temperature 
Tj> from equation ( jO] ). The density of states N?{E F ) at 
the Fermi level is determined exactly for the free-electron 
bands given in Figs. || and [|, the Debye temperature 9 
is arbitrarily put equal to 300-ftT and the effective spin- 
phonon interaction V is chosen in such a way that, in 
each case, the maximum transition temperature amounts 
to about 10K . A variation of V or 9 does not modify the 
positions of the maxima of T c noticeably. 

Changing the Fermi energy Ep stepwise over the entire 
energy range of the considered set of s-d bands, I calcu- 
late for each Ep first the total number n of electrons per 
atom by integration and then, for each i, the density of 
states Ni(Ep) of the ith band at the Fermi level. The cal- 
culated n is referred to as the average number of valence 
electrons per atom. 

The density of states NJ(Ep) is given by 



(3.1) 



Table [II gives the representations of a set of s-d bands 
in the hep structure as determined by the methods de- 
rived in Refs. [l^ and [jl]. Now the set consists of twelve 
bands since there are /i = 2 atoms in the unit cell. 

[T he representations given in T; 
Eq. (|b|) if we put d(a) = Xo(a)+X2(«), where Xo{ a ) 



III] satisfy 
1 

and X2(c>0 denote the characters of a in the representa- 
tions Dq and D2, respectively, of the three-dimensional 



where 



1 if the ith band belongs 

to the jth (T-band complex 
else. 



(3.2) 



In the following, we put j — U and j = L for the upper 
and lower maximum of T c , respectively. 
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FIG. 3. Free-electron s-d energy bands which serve as approximation of the energy bands of the hep solid solution alloys of 
the transition elements. The bands are numbered arbitrarily from 1 to 12 with the related points M of the reciprocal lattice 
being given in Table [jjb. Band 5 and 6 (dotted line) as well as band 7 and 8 (bold line) form superconducting band complexes 
of two bands each. Besides this, band 1 and 2 as well as band 3 and 4 form further superconducting band complexes not 
considered in this paper. 



A. Body-centered cubic structure 



For an integer number n of valence electrons, the. bec 
solid solution alloys studied by Hulm and Blaughero are 
elements with the space group G = and one atom per 
unit cell, cr-band complexes in this structure are simple 
a bands consisting of only one energy band. 

If n is not integer, the space group of the alloy is a sub- 
group G' of C\. An energy band recognized as a band in 
0\ stays a cr-band (complex) in the subgroup G', i.e., it 
stays a cr-band (complex) when the symmetry is lowered 
on the fixed atomic bec lattice. On the other hand, an 
energy band not forming a a band in 0% may become a 
cr-band complex when it is folded into the Brillouin zone 
of the subgroup G' . 



1. Lower maximum ofT c 



Table IV lists the representations of all the a bands 
in the bec structure. These representations are so-called 
extra or double-valued representations of the group Gj: 



of the wave vector k. 

The only a band which can be found among the energy 
bands in Fig. || is band 3 denoted by the dotted line. This 
cr band can be identified by means of a compatibility ta- 
ble for double-valued and single-valued representations. 
From the tables in Ref. |l8| we get 

D 1/2 x4 = r+ + r+ 

D 1/2 xH 25 =H++H+, 
D 1/2 xN 2 =N+, 

D 1/2 xP 4 = P 7 + P 8 . (3.3) 

These equations show that the s pin dependent Bloch 
functions 4>^ qm (r,t) [given in Eq. (C2)] of a band with 
the single- valued representations 

r' 25 , H 25 , Pi, and N 2 (3.4) 

can be transformed in such a way that at each symme- 
try point two functions form a basis of the double- valued 
representations 



r+ H+, P 7 , and N+ 



(3.5) 
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given in the second line of Table IV. Band 3 in Fig. || is 
a a band because the representations (3.4) can be found 
among the representations related to band 3. 
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FIG. 4. Transition temperature T c versus number n of va- 
lence electrons per atom in bcc solid sulution alloys of the 
transition elements, as calculated from Eq. ( [l.l| ) within the 
free-electron model by adjusting V such that the maximum 
T c amounts to about 10K. The positions of the maxima are 
independent of V. 

Lower maximum at n ~ 4.9 related to the superconducting 
band 3 denoted by the dotted line in Fig. 1 of part I and 
Fig. | 



Fig. |^ shows T^, as calculated with 
N£(E F ) = N 3 (E F ), 



(3.6) 



as function of the number n of valence electrons per atom. 
The maximum at n « 4.9 coincides with surprising ac- 
curacy with the experimental value n ~ 4.7 of Hulm and 
Blaugher.0 Obviously, the T c (n) curve is determined to a 
large extend by the symmetry structure of the s-d bands 
rather than by the real values of the energy E(k). This 
idea is corroborated by the calculated band structure of 
Nb [see Fig . 1 of part I], showing that a band with the 
symmetry ( |3.4| ) is roughly half-filled. This band has al- 
ready been considered in Fig. 1 of Ref. 



2. Upper maximum of T c 

The upper maximum of T c is not present in the bcc 
structure. Band 4, as denoted in Fig. |2]by the bold line, 
is characterized by the (double-valued) representations 



H+, P H 



and 7V 5 + , 



(3.7) 



indicating that at point P this band is connected with 
other bands, since the representation D1/2 x P 3 — P & is 
four-dimensional. Therefore, Pg, cannot be found in Ta- 
ble |v| and, hence, band 4 does not form a a band. This 
result, too, is confirmed by Hulm and Blaugher since the 



curve of the upper maximum does not go back to the in- 
teger value n — 6; see Fig. [|. Later measurements have 
shown that both Mo and W (having six valence electrons 
per atom) become superconducting at the relatively low 
temperatures 0.92K and 0.0012X, respectively. Their su- 
perconducting state is clearly related to band 3 because 
both metals have a band with the symmetry (3.4) in their 
band structure which is far from half, but not completely 
filled; see Fig. 1 in Ref. Ill 
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FIG. 5. As Fig. but the upper maximum at n m 6.5 
related to band 4 denoted by the bold line in Fig. ^. Be- 
cause band 4 is not a superconducting band within the bcc 
structure, this maximum does not exist in the bcc transition 
elements, while it exists in some ordered structures in bcc 
solid solution alloys. 

Fig. H shows a plot of T^(n) as calculated with 



+ 



+ 



N${E F ) = N 4 (E F ), 



(3.8) 



e.g., this curve is calculated by assuming that, neverthe- 
less, band 4 is a a band. The calculated maximum of 
T c (n) at n ~ 6.5 coincides with the experimental value 
n w 6.4 with high accuracy. This fact together with the 
steep rise of the experimental T c (n) curve for n > 6 in 
Fig. [I] suggests that band 4 becomes a cr-band complex 
when the ideal bcc symmetry of the electron system is 
disturbed, i.e., when band 4 is folded into the Brillouin 
zone of the space group G' of the Mo-Re or W-Re alloy. 

Unfortunately, in the literature I have not found any 
details of the symmetry of this alloy in the bcc solid solu- 
tion range. Possibly, the atoms are completely irregularly 
distributed over the bcc lattice sites. In this case, band 
4 becomes clearly a a band since the problem with the 
P3 state disappears. More likely, however, the atoms in 
these alloys are ordered or have at least a short-range 
order. Therefore, I have examined band 4 with the sym- 
metry (3.7) in four ordered structures in a body-centered 
cubic solid solution. 
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FIG. 6. Band 4 in Fig. g folded into the Brillouin zone 
of the CsCl structure, with symmetry notations as used in 
Table [v|. In this structure, band 4 does not form a supercon- 
ducting band complex (cf. Table |yj). 



First, consider the CsCl structure with the cubic space 
group 0\ and two atoms per unit cell. The double- valued 
representations of the four possible cr-band complexes (of 
two bands each) in the CsCl structure are listed in Ta- 
ble [v| Fig. H shows band 4 folded into the Brillouion 
zone of O x h . This set of two bands does not form a cr- 
band complex because the representations at point R, 
namely R$ , cannot be found in Table [v| The same re- 
sult we get for the NaTl structure with the space group 
7 h . Consequently, we can exclude that the CsCl as well 
as the NaTl structure is realized in bcc solid solutions of 
Mo-Re or W-Re since in both structures the measured 
upper maximum of T c cannot exist. 



The MoSi2 structure with the tetragonal space group 
D^ h has been observed rather frequently, but often with 
distortions of the atomic bcc lattice.Ell Table |Vl] lists the 
(double-valued) representations of the four possible a- 
band complexes in the MoSi2 structure and Fig. ^ shows 
band 4 folded into the Brillouin zone of DjjJ^. The re- 
sulting set of three bands forms a er-band complex as 
indicated by the second row of Table [v| 




( 



[P 6 ]P 7 
p 6 p 7 



n*,n; 




ri 



r 



r 



FIG. 7. Band 4 in Fig. || folded into the Brillouin zone 
of the MoSi2 structure, with symmetry notations as used in 
Table |vi]. According to Table yJ\, band 4 forms a supercon- 
ducting band complex in this structure. 

To summarize, in all the considered structures on an 
atomic bcc lattice, except for the CsCl and NaTl struc- 
ture, the upper maximum of the transition temperature 
as depicted in Fig. ^ may exist. 



B. Hexagonal close- packed structure 



In the hep structure we have two atoms in the unit 
cell and, hence, any cr-band complex consists of two en- 
ergy bands. All the possibl e cr- band complexes in this 
structure are listed in Table VII . By means of the com- 



patibility tables in Ref. 18 it can be verified that the set 
of band 5 and 6 as well as the set of band 7 and 8 de- 
noted in Fig. H by the dotted and bold lines, respectively, 
form cr-band complexes .t3 Figs. || and || show T r L an d 



respectively, as function of n calculated by Eq. (1.1) with 



NZ(E F ) = N 5 (E F ) + N 6 (E F ) 



(3.9) 



and 



Other conceivable structures on an atomic bcc lattice 
have the hexagonal space group D^ d , as, e.g., the Cdl2 
structure. [The International Tables give ±(-|,|,z) for 
the positions of the I atoms. The atomic lattice has 
bcc structure for z = | .] A hexagonal space grovui is 
suggested by the neighbouring "pseudo hexagonal" Eil a 
structure (see Fig. Q). It turns out that, as in the MoSi 2 
structure, band 4 becomes a cr-band complex in the Cdl 2 
structure when it is folded into the Brillouin zone of D^ d . 



N£{E F ) = N 7 {E F ) + N S (E F ). 



(3.10) 



Again the two maxima of T c at n ss 4.7 and n w 6.7 coin- 
cide almost exactly with the measured values of n ~ 4.7 
and n « 6.4, respectively. It should be noted that other 
possible definitions of the single free-electron bands yield 
slightly shifted positions of the max ima. This paper de- 
fines the bands as given in Eq. (2.1). 
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FIG. 8. Transition temperature T c versus number n of va- 
lence electrons per atom in hep solid sulution alloys of the 
transition elements, as calculated from Eq. within the 

free-electron model by adjusting V such that the maximum 
T c amounts about to 10K. The positions of the maxima are 
independent of V. 

Lower maximum at n ~ 4.7 related to the superconducting 
band complex consisting of band 5 and 6 which in Fig. ^ is 
denoted by the dotted line. 



IV. DISCUSSION 



Within the nonadiabatic Heisenberg model, the elec- 
trons of a narrow cr-band complex are constrained to form 
Copper pairs at zero temperature because in any normal 
conducting state the conservation of crystal-spin angular 
momentum would be violated; see part I of this paper. 
On the basis of this modified BCS mechanism of Cooper 
pair formation, a theoretical interpretation of the precise 
measurements of Hulm and BlaugherB, see Fig. [I], be- 
comes possible within an appropriate free-electron model 
and within the BCS theory. 

Hence, there is great evidence that the BCS theory 
yields the absolute energy minimum in the Hilbert space 
only if the formation of Cooper pairs is additionally con- 
strained by the constraining forces existing in a narrow, 
partly filled cr-band complex. Under this assumption, 
relevant results (but, of course, not all the results) of 
the strong coupling theory of superconductivity can al- 
ready be obtained within the weak coupling (BCS) limit. 
Within the approximations of this paper, I was able to 
calculate the numbers n of valence electrons per atom 
related to the maxima of the transition temperature T c , 
but not the values of T c . 

The calculated positions of the maxima of T c at va- 
lence numbers approximately equal to 4.9 and 6.5 in bec 
solid solution (see Figs. ^ and ||) are in good agreement 
with the measured values of 4.7 and 6.4 by Hulm and 
Blaugher; see Fig. [j] Obviously, the positions of the 
maxima of T c are determined to a large extend by the 
symmetry structure of the s-d electrons rather than by 
the special form of the periodic potential of the ion cores. 



This theoretical result emphasizes the universal character 
of the Matthias rule. 

In order to test the influence of the atomic structure 
on the transition temperature, the related maxima have 
been calculated in hep solid solutions, too. Here we 
get nearly the same positions as in bec solid solutions, 
namely n ~ 4.7 and n w 6.7; see Figs. || and |[ This re- 
sult suggests that the Matthias rule reflects a property of 
s-d electrons which is not influenced much by the lattice 
symmetry. 
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FIG. 9. As Fig. |S| but the upper maximum at n m 6.7 re- 
lated to the superconducting band complex consisting of band 
7 and 8 which in Fig. H is denoted by the bold line. 



This statement, however, should be modified: A too 
high symmetry of the lattice, i.e., of the periodic po- 
tential, may suppress superconductivity. A remarkable 
finding is that the upper maximum of T c at n « 6.5 ex- 
ists neither in the transition elements with bec structure 
nor in alloys with CsCl or NaTl structure, while it may 
occur in several ordered structures in a bec solid solution 
alloy of two components. These are, e.g., the tetragonal 
MoSi2 structure and several atomic arrangements with 
the hexagonal space group D\ h . Obviously, this fact is 
responsible for the abrupt increase of the T c curve in 
Fig. [I] between n = 6 and n = 6.3. 

The upper maximum in Fig. [l] is markedly narrower 
than the lower. This might indicate that a (short-range) 
order with the composition MoRe2 with n rs 6.33 is 
present in the solid solution. I believe that the ordered 
structure has the highest T c , though, at this stage, I can- 
not give a theoretical foundation for this speculation. 

The actual band structures of the transition metals 
contain also p-like Bloch states. In this paper I have 
taken into account p-like states in the bee structure but 
have ignored them in the hep structure; see Sec. II B. 
Therefore, I cannot exclude that the calculated maxima 
of T c in the hep structure are considerably shifted in re- 
ality. At this stage, I only can state that a pure s-d 
band structure in hep solid solutions leads to positions 
of the maxima which are similar to what is found in bec 
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materials. 



APPENDIX B: SUPERCONDUCTING BAND 
COMPLEXES 
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APPENDIX A: GENERAL 

Consider a metal with the space group G, the point 
group Go, and \i atoms in the unit cell. The positions of 
the atoms are written as 



T = R+g ll 



(Al) 



where R and Qi (i = 1 to /i) denote the lattice points 
and the positions of the ith atom within the unit cell, 
respectively. The elements 



= {«!*> 



(A2) 



of G consist of a point group operation a and a transla- 
tion 



t = f(a) + R, 



(A3) 



which is the sum of the nonprimitivc translation f(a) 
associated with a and a primitive translation R. 

The operators P(a) act on a function of position, /(r), 
according to 

P(a)f(r) = f(a- 1 f-a- 1 ?), (A4) 

and on Pauli's spin functions, 

u s {t) = 8 sU (A5) 

according to 

P(a)u s (t) = J2 d s's(a)u s >(t). (A6) 

The matrices [<i s < s (a)] are the representatives of the two- 
dimensional double- valued representation P1/2 of the 
three-dimensional rotation group 0(3), and s = ±^ and 
t = ±i denote the spin quantum number and the spin 
coordinate, respectively. _ 
The effect of K is given by the equations!^ 



and 



witl 



Ku s (t) = g g u- s (t), 
9±i/2 = T-i- 



(A7) 
(A8) 
(A9) 



The Bloch functions belonging to a cr-band complex 
can be unitarily transformed into spin dependent Wfs 
which are best localized, symmetry adapted to the para- 
magnetic group 



M 



G + KG, 



(Bl) 



and situated at the atoms in the sence that there is at 
each atom the center of symmetry of exactly one Wan- 
nier function. G denotes the space group and K stands 
for the operator of time inversion. A cr-band complex 
contains just as many bands as there are atoms in the 
unit cell. 

In this section, I shall give the group-theoretical con- 
dition to identify such a cr-band complex in a given cal- 
culated band structure. 

Consider a set of [i bands in the given metal and as- 
sume the Bloch functions in the points of symmetry P^ 
to be labeled by double-valued 2/x-dimensional represen- 
tations Ri of the group Gj of Pg. Ri may be irreducible 
or is the direct sum of irreducible double-valued repre- 
sentations. This set of fi bands forms a cr-band complex 
if two conditions are fulfiled. 

(1) In a cr-band complex, only representations Ri occur 
which can be written as a Kronecker product of a single- 
valued representation R^ with the representation D 1 / 2 , 



Ri = D 1/2 x R- k . 



(B2) 



Again, i?£ may be irreducible or is the direct sum of ir- 
reducible single- valued representations. 
(2) In a cr-band complex, the matrix rep resentatives 
D%(a) of the representation R^ in Eq. (B2) satisfy, at 
each symmetry point Pg, the equation 

trace Dj:(a) 

= d(a)e- la ^ r J2 riii^e-^- 01 ^ for a £ G % , (B3) 



where d(a) stands for the representatives of any real one- 
dimensional single- valued representation, say To, of Go 
and 



rii(a) 



1 if (XQi + t = Qi 

else 



R 



(B4) 



with R being a primitive translation (which, of course, 
may be equal zero). 

Eq. (fH) is derived from Eqs. (1.8) and (4.28) of Ref. |l| 
(cf. Sec. 3 of Ref. [l6|) assuming that there is one Wannier 
function with To symmetry at each atom. Eq. ( |B3| ) is sat- 
isfied for all a £ Gg if it is satisfied for one representant 
of each class of G^. 
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APPENDIX C: SPIN DEPENDENT WANNIER 
FUNCTIONS 

In this section the spin dependent Wfs belonging to er- 
band complexes will be defined. Consider a cr-band com- 
plex consisting of fi bands. The Bloch functions <fj: q (r) 
of this band complex can be unitarily transformed into 
spin dependent Wfs 

w im (r- R- Qi,t) 

BZ M 

where the spin dependent Bloch functions 4>j: qm ('r, t) have 
/c-dependent spin directions, 



kqra 



(r,t)= fsni(q,k)u s (t)ip nq (r). 



(C2) 



The first sum in Eq. (|Cl|) runs over the N wave vectors 
of the Brillouin zone (BZ), q = 1 to /i denotes the band 
index, the /x-dimensional matrices [gi q {k)\ are, for each 
fc, unitary, and the two-dimensional matrices [f S m{q^k)\ 
are, for each q and fc, unitary, i = 1 to fx labels the 
Wfs at different atoms in the unit cell and the crystal- 
spin quantum number m — ±i distinguishes between the 
functions at the same atom. 



In a cr-band complex, that means, if Eqs. (B2) and 

a q (k)} and [/ s . 



(B3) are satisfied, the matrices. 1$ 
can be chosen in such a wayEffijEfl that 

(1) the spin dependent Bloch functions ^ > g gm (^')*) vary 

smoothly through the whole k space and 

(2) the Wfs are symmetry adapted to M p according to 

= di(a) ^ d m ' m (a)wf, m ,(r,t) for a € G (C3) 



and 

K Wf m (r,t) = ±w f _ m (r,t), 
with the abbreviations 



i(r-R - Qi,t) 



and 



f ' = a(R + g t ) + t. 



(C4) 



(C5) 



(C6) 



\ and 



The plus in Eq. (C4) is defined to belong to m = 
the minus to m = — \. As a consequence of Eq. (B2), the 
two-dimensional matrices [d m > m (a)] belong to the repre- 
sentation -D1/2 of 0(3) and the coefficients di(a) are c- 
numbers with |cZ^ (o) | = 1. The smoothness of the spin 



dependent Bloch functions ^>y. qm if-, t) guarantees that the 
spin dependent Wfs are optimally localizable. 

1, the c- numbers di(a) are 



[If in Eq. (|BJ) m{a) 



equal to d(a) as given in Eq. (B3). In all cases we have 
di(a) — Dji(a) with j labeling the non- vanishing ele- 
ment in the ith column of the matrix [Dji(a)] given in 
Eq. (2.18) of Rcf. |l6| which has only one non-vanishing 
element in each column.] 
If in Eq. ( |C2| ) we have 



fsm(q, k) = S s 



(C7) 



the Wfs are usual Wfs. We have spin dependent Wfs if 
the matrix [f sm (q, k)] is not independent of k. 

Usual Wfs could be constructed only if the cr-band 
complex would not be connected with other bands (if 
we still demand Lhat. the Wfs are best localized and sym- 
metry adapted). In the metals, however, any energy 
band is connected which other bands at the point and 
lines of symmetry. Therefore, the matrix [f S m(q, k)] can- 
not be chosen independent of k. 

The abbreviation in Eq. ( |C5| ) is possible since there is 
one Wannier function at each atom in a cr-band complex. 
Therefore, the description of the nonadiabatic Heisenberg 
model given in Sec. II of part I of this paper need not 
be changed in the case of cr-band complexes consisting of 
more than one band. 
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(a) 



(b) 



1 ( 0, 0, 0) 

2 ( 0, 0, 1) 

3 ( 1, 0, 0) 



4 (0,- 

5 (0, 

6 (0,- 



1) 
-1) 
0) 



1 ( 0, 0, 0) 

2 ( 0, 0, 1) 

3 ( 0, 1, 0) 

4 ( 0, 1, 1) 

5 (-1, 1, 0) 

6 (-1, 1, 1) 



7 ( 1, 0, 0) 

8 ( 1, 0, 1) 

9 (-1,0,0) 

10 (-1, 0, 1) 

11 ( 1,-1,0) 

12 ( 1,-1, 1) 



TABLE I. The points M of the reciprocal lattice related, 
according to Eq. ( |2.l|) , to the ith free-electron energy band 
given in Fig.Ji] and |5[ respectively. Each column lists i to- 
gether with Km in the basis given in Table 3.3 of Ref. [l]. 
(a) bcc structure, (b) hep structure. 



2H 



2r- + r+ + r- + r+ + r 6 - 



2H X + H 2 + 3H 3 
K 2 + K 3 + Ki + 3K 5 + K 6 
2A 1 + 2A 3 
4Li + 2L 2 



3M+ + M{ + 



3M; 



Mi 



Ml 



M, 



•AT; 



TABLE III. Representations of an ideal set of s-d bands in 
the hep structure [as determined by means of Eq. (|B3|), see 
Sec. 1IB | with notations given in Sec. A3-1 of Ref. 



n 



ri + ri2 

Hi + H12 

P1+P3 + Pi 
ZNi + N 2 + N 3 + N A 



25 
#25 



TABLE II. Representations of an ideal set of s-d bands in 
the bcc structure with notations given in Sec. A3-3 of Ref. [l^. 



r+ 




Ps 


K 




H+ 


Pt 


K 


r 6 " 


Hs 


Pe 




r 7 " 


H 7 - 


Pi 


N 5 ~ 



TABLE IV. Double- valued representations R? of the four 
superconducting bands in the bcc structure with notations of 
Ref. |l| 



2r+ 


^6 


+ Rj 


M+- 






fX 6 " 


2Tf 


R+ 


+ Re 


M+- 






f Xj 


2r 6 ~ 


i?+ 


+ Re 


M 6 "- 


h Mf 




f A 6 - 


217 


Rt 


+ Rj 


M 6 "- 


l-Mf 


x+- 


f Xj 



TABLE V. Double-valued representations Ri of the four superconducting band complexes (of two bands each) in the CsCl 
structure with notations of Ref. hq. The CsCl structure has the cubic space group 0\. 



2r+ 


+ r 6 " 


2Z+ 


+ Z 6 " 


2X+ 


+ * 5 ~ 


2P 6 


+ P7 


2N 3 


+ 2N+ 


+ N 3~ 


+ iv 4 ~ 


2r^ 


+ r-7 


2Z+ 


+ Zj 


2X+ 


+ *s 


P 6 ^ 


2P 7 


2N+ 


+ 2N+ 


+ iv 3 - 




r«T + 


217 




-2Z 6 " 




-2Xr 


P 6 4 


2P 7 


N+ 4 


-N+ + 


2N s ^ 


-2iv 4 - 


r+4 


217 




-2Zf 




-2X 5 - 


2P e 


+ P7 


N+-\ 


-N+ + 


2^ ^ 


~2N± 



TABLE VI. Double-valued representations Ri of the four superconducting band complexes (of three bands each) in the 
MoSi2 structure with notations given in Table 6.14 of Ref. Il3l The MoSi2 structure has the tetragonal space group D^ h . 
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Tf+Ta H 5 + H- + H g K 8 + K 9 A e L 3 + L 4 M+ + M 5 " 
Tf+r+ H 4 + H 6 + H s K? + K 9 A e L 3 + L± M+ + M 5 ~ 

TABLE VII. Double-valued representations Ri of the two superconducting band complexes (of two bands each) in the hep 
structure with notations given in Ref. \m. 
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